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We propose an optoelectronic scheme to define and manipulate an indirect neutral exciton qubit 
within a quantum dot molecule. We demonstrate coherent dynamics of indirect excitons resilient 
against decoherence effects, including direct exciton spontaneous recombination. For molecules with 
large interdot separation, the exciton dressed spectrum yields an often overlooked avoided crossing 
between spatially indirect exciton states. Effective two level system Hamiltonians are extracted by 
Feshbach projection over the multilevel exciton configurations. An adiabatic manipulation of the 
qubit states is devised using time dependent electric field sweeps. The exciton dynamics yields the 
necessary conditions for qubit initialization and near unitary rotations in the picosecond time scale, 
driven by the system internal dynamics. Despite the strong influence of laser excitation, charge 
tunneling, and interdot dipole-dipole interactions, the effective relaxation time of indirect excitons 
is much longer than the direct exciton spontaneous recombination time, rendering indirect excitons 
as potential elemental qubits in more complex schemes. 

PACS numbers: 71.35. Gg,71.35.-y, 03.67.-a, 73. 21. La, 78.67.Hc 
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I. INTRODUCTION 

Semiconductor quantum dot molecules (QDMs) are 
potential building blocks for solid state quantum com- 
putation architectures!^ These devices have remarkable 
electronic and optical pr oper ties, arising from their 
molecular exciton spectrumPEl Moreover, the optical re- 
sponse of QDMs is highly tunable using external elec- 
tric and magnetic fields. Quantum optics techniques 
have demonstrated coherent phenomena in QDMs, such 
as Rabi oscillations and level anticrossing of excitonic 
dressed states.^HS Strong localization of charge and spin 
in these structures permits multiple ways of harness- 
ing charge and spins qubits, all limited by uncontrolled 
interactions within the molecule environment P On one 
hand, spin qubits in self-assembled quantum dot struc- 
tures have received a great deal of attention due to 
the large coherence time of spins localized in quantum 
dots, limited mostly by their hyperfine interactions with 
the QDM nuclear spin reservoirP In particular, opti- 
cal spin initialization and non-destructive measurements 
have already been implemented in QDMs. 10 On the other 
hand, charge qubits have typically shorter decoherence 
times, limited by spontaneous exciton recombination and 
electron-phonon interactions!^ The latter limitation can 
be largely suppressed at very low temperatures while the 
former is more subtle. For neutral spatially direct exciton 
qubits, with logic states typically embodied in the pres- 
ence or absence of an exciton in a single QD, spontaneous 
recombination is fast (~ Ins) and highly detrimental, due 
to the large direct exciton oscillator strength. 

In this work we investigate the exciton dynamics 
of optically driven and electrically gated QDMs cou- 
pled by c harg e tunneling and Forster energy transfer 
(FRET)PHHT We argue, theoretically and numerically, 
that an exciton dressed qubit EuGZQ w ith logical states 



constituted by neutral indirect excitons, can be effec- 
tively extracted from the QDM exciton dressed spectrum. 
Our work indicates that a control scheme is devisable us- 
ing external electric field sweeps. By these means, the 
qubit can be initialized and rotated multiple times with 
high fidelity, well beyond the spontaneous exciton recom- 
bination time scale. Furthermore, we devise a read out 
scheme, using an adiabatic population transfer of the out- 
put indirect state into an auxiliary direct exciton state. 
This opens the possibility of a realistic realization of a 
neutral exciton qubit with enhanced characteristics and 
subsequent coherent manipulation using optical and elec- 
trical means. For typical QDM structures, we find that 
one can anticipate relaxation/decoherence times of at 
least two-orders of magnitude larger than spontaneous 
recombination times. Further separating the dots in the 
QDM could give rise to even longer coherence times, mak- 
ing them suitable for implementing complex multiqubit 
architectures. 

In Sec. [IT] we introduce a realistic model for the QDM, 
that takes into account all relevant electron-hole states, 
and processes at the relevant energies. An excitonic 
dressed spectrum and population bias map is introduced 
in Sec. IIIII and used to indicate the different molecular 



resonances (level anticrossings). In Sec. IV we employ 
a Feshbach projection formalism, 18 adiabatically elimi- 
nating selected exciton transitions and extracting an ef- 
fective Hamiltonian describing locally the relevant level 
anticrossings. In particular, we demonstrate that adia- 
batic elimination of the direct exciton transitions leads 
to an effective qubit subspace consisting of two long lived 
indirect excitons. In Sec. |IV A] we describe the qubit dy- 
namics in two regimes, strong vacuum-indirect exciton 
coupling (qubit initialization) and strong coupling among 
two molecular indirect excitons (qubit rotation) . We dis- 
cuss the role of FRET and biexciton states in Sec. IIVBI 



and find that the qubit subspace is resilient against their 
detrimental effects as well as to small corrections to the 
laser detuning. Sec. [V] introduces an applied adiabatic 
bias ramp that implements the initialization, rotation 
and readout of the qubit logical states. It is demonstrated 
that the qubit can be initialized and rotated with near 
unity fidelity within a picosecond scale. Finally, in Sec. 
VTJ we develop a method to extract effective decay rates 
of the molecular indirect excitons during initialization 
and rotation regimes, which results in relaxation times 
many orders of magnitude larger than the manipulation 
times. 



II. MODEL 

The QDM consists of two vertically stacked non- 
identical "top" (T) and "bottom" (B) quantum dots. 
The dots are separated by a barrier of thickness d and 
subject to an applied axial electric field F that results 
from the application of a top gate voltage. This is re- 
alized by placing the QDM in a n-i Schottky junction. 2 
The QDM is pumped by a broad square laser pulse of fre- 
quency w, which may excite different nearby exciton lev- 
els. The pulse duration is long enough, typically > lps, 
to capture several amplitude oscillations of the excitonic 
populations. 
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FIG. 1. (Color online) (a) Bare exciton level configuration 
corresponding to Hamiltonian Eq. (TTj) . Solid arrows indicate 
exciton transitions mediated by different processes of tunnel- 
ing (t e ,th), optical pumping (fi) and Forster transfer (Vf); 
dashed arrows represent radiative channels, (b) Effective level 
configuration according to Eq. [7] after adiabatic elimination 
of the direct transitions; the tilde on labels represent effective 
couplings and decay rates. 

We denote exciton bare states by \ e h Be h T X) 7 where 
e_B(T):^s(T) = {0,1} are the electron and hole occu- 
pation numbers, resulting in a total of five states, as 
shown in Fig. [T^,. The basis of this excitonic Hilbert 
space contains: the vacuum |qo-^0; two single direct exci- 
ton states, \\qX) (bottom exciton) and loi-X") (top ex- 
citon); and two single spatially indirect exciton states 
loi-X") and \%X). The Hamiltonian in the rotating wave 



approximatiorPSfli] j s given by 
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where the columns are associated with the states IoqX), 
\%{X), \f x X), \%X), and \{%X). Diagonal matrix el- 
ements represent detunings of the exciton levels from 
the laser energy. An applied axial electric field F re- 
sults in Stark shifts for the indirect exciton detunings 
<5io(oi) = eiojf(oix) - ftw given by 
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(2) 



Likewise, the excitation pulse electric field envelope gen- 
erates an optical matrix element, fi = (qqX\jX-e\ 1( S 01 Ix), 

for direct exciton transitions! 2 ^ 2 ^ Matrix elements, t e ,th, 
describe single particle interdot tunnelings for electron 
and hole respectively.^ Vf is an interdot "hopping" for 
an entire electron-hole pair, arising from the dipole-dipole 
interaction coupling two direct exciton states, 



V F = 



4ireoe r d 3 



(3) 



where e r is the dielectric constant, and ^t(b) ~ 6.2eA are 
the interband transition dipole moments. These are as- 
sumed parallel to each other (and perpendicular to their 
separation), which gives an orientation factor k ~ l! 2 ^ 
In the numerical calculations we use d ~ 8.4nm, which 
yields a value of Vf = 80/ieV. This interaction allows 
the Forster energy transfer mechanism by which a donor 
QD transfers its exciton energy to the acceptor neighbor- 
ing dot, effectively resulting in the non-radiative interdot 
"hoping" of the excitonP^lE] 



We also consider radiative decay rates, T 



x 



TX 



Ins, describing the spontaneous recombination of exci- 
tons with spatially direct characterPS To that effect, the 
exciton dynamics is obtained from solutions to the Lind- 
blad master equation, which yields the time evolution for 
the density matrix of the systempSHZD 
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L(p). 



(4) 



The first term on the right describes the coherent evo- 
lution of the excitation dynamics, H(t) being the full 
Hamiltonian of the system, Eq. (fTl). The second term, 
L(p), incorporates dissipation processes, 



L(P) 
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where Pj = \j)(j\, and \j) is an exciton which relaxes 
into a state \i), with rate T^-. The dynamics requires 
the solution of iV 2 coupled differential equations for a 
./V-dimensional Hilbert space. 



As we will explain in detail later, additional states con- 
sidered in the model correspond to neutral biexcitonic 
states, Sec. |IVB| whereas exciton states arising from ex- 
cited states of the electron (hole) and charged excitons 
are not considered (assumed to be far removed from the 
manifold of interest) . We also assume that charge tunnel- 
ing into the contact reservoirs and spin-orbit interactions 
are negligible, rendering electron-hole exchange decoher- 
ence processes unimportant, see Sec. |VD| 



III. EXCITONIC DRESSED SPECTRUM 

The interplay of charge tunneling, incident radiation 
field, and Coulomb interactions, results in coherent in- 
terdot coupling. This coupling yields complex molecu- 
lar states that are superpositions known as the dressed 
excitonic states. Only a subset of these superpositions 
leads to allowed transitions, with field dependent ampli- 
tudes and energies, which results in anticrossings in the 
dressed exciton spectra of the QDM. The time evolu- 
tion of molecular states under optical pumping results in 
Rabi oscillations, which can be time-integrated to yield 
the average occupation of the excitons involved in the 
molecular states. The occupation of such excitons can 
be probed in principle by differential transmission of a 
weak probe measuring the population of a particular ex- 
citon. It is then possible to construct a level anticrossing 
(LACS) population map, using the time integrated dy- 
namics of the Hamiltonian, Eq. (nl. For a given map 
coordinate, (F,ftw), the integrated population is given 

by Pi = (l/*i) J L Pu{t)dt, where tz, ~ 500ps stands 
for the constant-amplitude pulse duration, long enough 
to capture several amplitude oscillations of the exciton 
populations; only a few Rabi oscillations are necessary 
to reliably compute pi. Then, any exciton state pop- 
ulated under pumping will exhibit a relative amplitude 
Pi(F, Hlo) and contribute to features on the corresponding 
map. At each coordinate, (F,hu), two or more excitons 
share population if they have non- vanishing components 
in the dressed state; then by examination of maps cor- 
responding to individual excitons, one can reconstruct 
the entire dressed spectrum of the system. Alternatively, 
one can compute the population map of the vacuum state 
\qqX = |0), so that the complete dressed LACS spectrum 
will correspond to all (F, huj) coordinates where this es- 
tate is depopulated; such is the case of Fig. 2. 

In contrast, Fig. 3a shows the QDM eigenvalue spec- 
trum as function of applied electric field F for a con- 
stant value of the pump laser energy, hu — Ei aser = 
1299. 6meV and system parameters as in Fig. 2. The top 
right box indicates a level anticrossing that has a corre- 
spondent optical signature in the population map of the 
vacuum state (top right box in Fig. 2). Anticrossings 
at zero energy between the vacuum IqqX) and exciton 
states loi-X") and \%X) appear at electric field values, 
F ~ —38.6 and 43.4kV/cm, respectively. These values 
coincide with the corresponding indirect exciton popu- 
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FIG. 2. (Color online) Level anticrossing population map of 
the QDM vacuum state |oo^}- As function of applied bias and 
pump laser energy the vacuum depopulates to other excitons 
at each resonance. Upper right dashed box indicates resonant 
excitation into an indirect exciton, %X, at Fi = 43.4kV/cm, 
for a laser energy Ei aS er = 1299. 6meV. The central anticross- 
ing mixes the two indirect excitons at Fr — 2.3kV/cm. See 
system parameters in Ref. |27] . 



lation signatures in upper part of Fig. 2, indicated by 
the dashed horizontal line. In a similar way the cen- 
tral box in Fig. 3a encircles a very narrow anticrossing 
between indirect excitons |oi-X") and lifjA), occurring at 
F ~ 2.3kV/cm. Notice that each anticrossing occurs 
energetically far away from each other, indicating that 
their eigenstates superpositions are only weakly coupled 
to the others. This suggest that an effective Hamiltonian, 
represented in the basis of the qubit subspace, should re- 
produce these anticrossing signatures. 



IV. QUBIT EXTRACTION 

In the spectroscopy of QDMs, the appearance of a level 
anticrossing signature points to the onset of an important 
interdot interaction. These interactions are in some cases 
not straightforwardly explained by the off diagonal ma- 
trix elements of the Hamiltonian, Hy, connecting two 
allowed states. For example, two excitons |i) and |j), 
with a very weak oscillator strength can couple strongly 
to the radiation field via higher order transitions medi- 
ated by non-optical processes, such as charge tunneling. 
Such is the case for indirect excitons in our model, which 
couple via higher order processes to the radiation field, 
yet they are assumed to have zero oscillator strength. 



These states have well-defined optical signatures in the 
LACS map in Fig. 2 and exhibit an anticrossing with the 
vacuum in the eigenvalue spectrum in Figs. 3b and 3d. 

This coupling of two indirect excitons via higher order 
transitions, and their respective coupling to the vacuum, 
should be revealed by an effective Hamiltonian that de- 
scribes the same physics as the original, but constrained 
to a sector of the Hilbert space whose wave functions 
correspond just to the eigenvalue spectrum in the anti- 
crossing region. This projected Hamiltonian is of reduced 
dimensionality and should have non-zero off-diagonal ma- 
trix elements connecting the states involved. We employ 
the Feshbach projection operator formalism, 18 which per- 
mits the derivation of an effective Hamiltonian with ex- 
citon eigenstates of pure indirect character. The direct 
exciton sector of the Hamiltonian is adiabatically elim- 
inated ("projected out"), and its dynamical effects be- 
come embedded in the matrix elements of the effective 
Hamiltonian. The requirements of adiabatic elimination 
are satisfied by two conditions: (1) adiabatic variation 
of all external fields (2) ability to isolate spectrally, by 
tuning the excitation energy (off-resonant condition) and 
applied electric field, the confluence of two excitons from 
the remaining exciton manifolds. The Hamiltonian ob- 
tained by projection, would describe the time evolution 
of qubits in the unitary and dissipative regimes. 

In the following discussion, for simplicity we consider 
a closed quantum system with the Hamiltonian given by 
Eq.fTJ The Hamiltonian can be separated in to two parts, 
H = Hq + V, where Hq is the unperturbed diagonal part, 
and V a perturbation. Let V be the relevant subspace 
spanned by the excitons that self-avoid at a chosen sys- 
tem resonance. In the same way let P and Q — 1 — P be 
projector operators onto and outside of V, respectively. 
The effective Hamiltonian is given bjff^ 



H(z) = PH P + PR(z)P , 



(6) 



with z — E±ie, where E and e are the real and imaginary 
parts of the complex energy eigenvalue z. The first term 
of H is the leading unperturbed part of the Hamiltonian 
inside V, with the second term containing the level shift 
operator, R(z) = V + VQ[z - QH a Q - QVQ^V, pro- 
jected onto V . The latter term can be seen as a Hamilto- 
nian that permits the calculation of the energy level shifts 
with respect to the unperturbed levels. Allowing the 
Hamiltonian to depend on its eigenvalues z, makes the 
eigenvalue equation non-linear. Additionally, analytic 
continuation of the eigenvalues into the complex plane 
allows the introduction of non-Hermitian Hamiltonians 
that incorporate dissipation processes taking place out- 
side the relevant subspace, V . Self consistent solutions 
to the non-linear eigenvalue equation are used to obtain 
the eigenvalue spectrum in the vicinity of a level crossing 
and anticrossing. Near a level anticrossing (and in the 
absence of accidental degeneracies) there is a unique self- 
consisting solution of z(F) for each value of the applied 
electric field F. 
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FIG. 3. (Color online) (a) Dressed exciton eigenvalue spec- 
trum at fixed excitation energy fko = 1299. 6meV. Boxes show 
anticrossings associated to couplings to the radiation field and 
between indirect excitons. These features have a one to one 
correspondence to the optical signatures in Fig. [2] (b) and 
(d) zoom into anticrossings near E — 0. (c) is the central an- 
ticrossing that mixes the indirect excitons spanning the qubit 
subspace. 



A. Indirect exciton qubit Hamiltonians 

In order to obtain the dynamics of molecular exci- 
tons with spatially indirect character, we project the 
total Hamiltonian (fil) onto the three level system, 
|°°A),|i°A),|°jA), shown in Fig. lb. The resulting 
Hamiltonian incorporates an effective coupling among 
the indirect excitons, U, and effective couplings of in- 
direct excitons to the radiation field, Ooi and fiio, re- 
spectively. We assume that the direct exciton levels 
are resonantly coupled by FRET, <5io = <5oi = A, and 
that the bottom and top QDs couple with the same 
strength to the radiation field, Q. This assumption 
yields, Ooi = fho = flj. Let P\ be the subspace sub- 
tended by the vacuum and the indirect exciton levels; 
then a projection of Hamiltonian (fTl) onto Pa gives 

Aoo(z) fli(z) &i(z) 

m A \z)= \ tii(z) A h a(z) + A s U(z) 

JM*) U(z) Aoi(z)-A s 

(7) 

The matrix elements of H A (z) contain level shift detun- 
ings Aj and effective couplings U and Qj. The indirect 
exciton effective detunings are given by 



Aio/on(z) = (Jiocoii + 5r(z) . 



(8) 



5j(z) 



(z-A)(t 2 e +tl) + 2t e t h V F 



(z-AY 



V 2 



(9) 



while the shift of the zero of energy (we have set Sqq = 0) 
is given by 



nn v / 



2n 2 



A + V F 



(10) 



correspondingly, the indirect exciton effective coupling to 
the radiation field arc both given by 



iM*) 



Ofte + t h ) 

A + Vf-z 



(11) 



Figures 3b and 3d show the anticrossing gaps opened by 
the couplings fli in Eq. (M). Both gaps have a width 
2f2j = 65/ieV, occurring at applied electric field values 
of F ~ -38.6 and 43.4kV/cm, respectively!^ 1 fi/ is di- 
rectly proportional to the direct transition dipole matrix 
element fi and tunneling amplitudes (t e + th). In other 
words, molecular indirect excitons are provided with an 
effective oscillator strength when the dots in the QDM 
are tunnel coupled. This leads to the possibility of reso- 
nant excitation of the indirect excitons, even if they have 
a vanishing intrinsic oscillator strength. This explains the 
"lighting up" of indirect excitons in Fig. 2 when optically 
driving the QDM at electric field values |.F| > 20kV/cm. 
The relevant interaction between two neutral indirect 
excitons is given by 



U(z) = 



2{z- 



A)t e t h + 
(z - A) 2 



(tl 



ti)v F 



V 2 
v F 



(12) 



which dominates for values of electric field Fd = 



2.3kV/cm. Eq. (12) shows that U is independent of the 
laser intensity embodied in the direct dipole matrix ele- 
ment fl. This represents an important desirable feature 
for a qubit defined in the indirect exciton subspace. It 
implies that for a fixed laser energy, the effective qubit 
subspace generated by the indirect excitons \oiX) and 
\%X) is effectively shielded against the external distur- 
bance of the intense optical field and less susceptible to 
the effects of spontaneous direct exciton recombination. 
On the other hand, U arises predominantly from elec- 
tron and hole tunneling, with a weak contribution from 
FRET; this means that the molecular indirect subspace 
evolves mainly by its internal dynamics. For the system 
under consideration, see Ref. [27], we find 2U ~ 45/ieV. 
A more concise qubit Hamiltonian is obtained by pro- 
jecting Eq. (l} onto Vi = {|o?X}, \%X)}. Then, one 
obtains a two-level Hamiltonian describing the spectrum 
at the central anticrossing when U dominates. The pro- 
jection results in 



m i >( z ) = 



HUz) U(z) 
U(z) H&(z) 



+ S(z)(<t x + I), (13) 



in terms of the matrix elements of Eq. ph , and correc- 
tions with ax and / being the x-Pauli and identity ma- 
trices, respectively. The correction term £(z) = z _£ , , 



describes the small perturbations arising from optical ex- 
citation and reflects the fact that the effective subspace 
is not perfectly isolated when the radiation field is on. 
Notice, however that the correction terms disappear in 
the absence of pumping (O = Oj = £ = 0). We also 
emphasize that H^\ defined on the subspace Pi, pro- 
vides a better description of the qubit rotation, while 
&( ' is more adequate and convenient for describing the 
initialization of the qubit system (via the effective cou- 
pling of indirect excitons to the light field). However, 
both subspaces and associated Hamiltonians are suitable 
to describe the intrinsic qubit dynamics. 



B. FRET and biexciton effects 

As mentioned above, when direct exciton transitions 
in the two dots are near resonant, the Forster energy 
transfer mechanism plays an impo rtant role on the di- 
rect exciton superradiant dynamics! 28 * 29 * Typically Vf — 
0.08meV, for interdot separation d ~ 8.4nm. This is a 
small value in comparison with electron tunneling £ e p31 
but appreciable enough to split the direct exciton spec- 
tral lines and redistribute the exciton population (spec- 
tral weight) among the molecular states in a steady state 
M 



regime. — Interestingly, the denominators in the effective 



Hamiltonian matrix elements in Eq. ( 11 ) and ( 12 ) exhibit 



a dependence on Vp as a correction to the direct exci- 
ton detuning A. Therefore, any influence of Vf can be 
strongly suppressed whenever A 3> Vf- In this regime, 
possible dephasing effects due to FRET would be sup- 
pressed as well. In our model, |A| ~ 51.5 meV, assur- 
ing the indirect exciton qubit subspace is indeed shielded 
against the perturbation effects of FRET. 

On the other hand, strong laser excitation can pump 
additional exciton levels outside the relevant subspace 
of consideration. The closest excitations are biexciton 
resonances, which cannot be in principle ignored in the 
dynamics of single excitons, as their detuning is at most 
a few meV.^l The pumping of biexcitons in either QD, 
||qX), lo^X), expands the bare exciton basis to 14 states 
\ e h Be h T X), with possible double occupancy of the single 
particle levels, es(T),/iB(T) = {0,1,2}, which becomes 
more significant with higher excitation power and/or 
short laser pulses. However, the detuning of the biex- 
citon levels, and the need for a direct exciton prior to its 
formation, result in weak perturbativc effects of the biex- 
citon level manifolds for the values of the matrix element 
n considered here. Moreover, the biexciton manifolds de- 
couple once the excitation power switches off during the 
dynamical control procedure, see Sec. [V] Other excita- 
tions, such as LO phonon resonances, appear ~ 35meV 
above the lowest exciton transition for GaAs, and can 
be safely ignored.^ We notice that the chosen structure 
parameters (QD confinement sizes) result in excited elec- 
tron and hole states (and associated excitons) far from 
the relevant anticrossing gaps, 2U and 2ilj, so that these 
other excitations can be safely ignored (see Sec. VD). 



V. COHERENT ROTATION OF INDIRECT 
EXCITON QUBIT 

In what follows, we consider the system as a fully open 
quantum system and consider explicitly radiative recom- 
bination of direct excitons. In this sense, we analyze our 
results in terms of numerical solutions to the Lindblad 
master equation, Eq. |4l with all 14 excitonic states in- 
cluded. Our discussion of the projected subspace Vi, 
indicates that the molecular indirect exciton subspace 
is indeed weakly influenced by interdot energy transfer 
mechanisms, Vp, and excitation power, Q. This sug- 
gests that we can achieve control of the indirect exci- 
ton qubits by tuning the effective coupling strengths Qi, 
U and by application of external time dependent elec- 
tric fields (we will discuss in detail dissipation processes 
below) P2 We use a cyclic adiabatic variation 



Sec. 



VI 



of the applied field, F(t), at fixed excitation energy Hu>, 
between a regime where the system effectively contains 
two-levels mixed by the coupling Qi, into a regime where 
the system contains two levels mixed by U. One can use 
short adiabatic bias pulses for qubit initialization and 
rotation operations.^ Figure 4a shows a cyclic sweep of 
applied bias, the left arm (0 < t < 0.17ns) indicates 
the initialization regime, shown in more detail in Fig. 4c. 
The slow forward bias ramp (0.17 < t < 1.22ns) drives 
the system into the qubit rotation regime, see Fig. 4d. 
The plateau in the bias pulse (1.22 < t < 2.99ns) corre- 
sponds to the rotation regime, and its tunable duration 
determines at which particular time one decides to ro- 
tate the input state (red curves) or not. If rotated, the 
reverse bias adiabatic ramp (2.99 < t < 4.05ns) transfers 
the output state population (blue curves) into a direct 
exciton (green curves) , which in turn depopulates subse- 
quently into the vacuum (dashed curves) for (t > 4.05)ns, 
see Fig. 4e. The overall dynamics of the implemented co- 
herent control is shown in Fig. 4b. In what follows, let 
us discuss each region in more detail. 



A. Initialization 

The couplings U and fli dominate in different field 
regimes, therefore the representation subspace of Hamil- 
tonian Eq. Q can be decoupled in three different regions. 
When 17/ dominates, for large values of electric field and 
positive energy detuning, we can construct two projected 
subspaces, spanned by the basis vectors {|qo)> loi-^Q} an d 
(loo)) I io-^O }> respectively. For an excitation energy of 
hio = 1299. 6meV, coherent Rabi oscillations are induced 
in each of these subspaces for applied electric field values 
of F = —38.6 and 43.4kV/cm, which corresponds to res- 
onant excitation at either level anticrossing, shown in fig- 
ures 3b and 3d, respectively. This allows the implemen- 
tation of 7r rotations within a time interval 



2f2r 



63.6ps, 

enabling the possibility of initializing the system in either 
of the logical states, \%X) or \%X), respectively. On 
the other hand, SI/ depends inversely on the direct exci- 



ton detuning, A, which sets an upper bound such that 
fli ^$> Tx , since otherwise the fidelity would be hampered 
by spontaneous recombination. 

Figure 4c shows the initialization of the indirect ex- 
citon \%X) by a 3n rotation. When driving the sys- 
tem at the anticrossing, corresponding to the coordinate 
(Fi = 43.4k V/cm, Hu = 1299.6meV), the initialization 
takes place after switching off the pulsed resonant excita- 
tion at a time t = 200ps. The initialization occurs with 
near unity fidelity, F = Q^X\p(U)\HX) ~ 0.97, due to 
the almost perfect isolation of the subspace Va- 
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FIG. 4. (Color online) Exciton population dynamics sub- 
ject to an applied cyclic bias pulse, (a) Applied bias pulse 
sweep of duration r ~ 5.5ns, the cycle sweeps the interval 
F = [43.4, 2.3]kV/cm. (b) Exciton population dynamics cor- 
responding to the cyclic sweep above. Notice (a) and (b) 
panels have the same time scale, while (c)-(e) show details of 
different regions, as indicated, (c) Rabi oscillations between 
looA") and |o?X) excitons. After a 37r rotation, at t — 200ps, 
the laser is turned off and qubit is initialized with near unity 
fidelity, (d) Rabi oscillations inside the qubit subspace. (e) 
Read out scheme of the output qubit state, via tunneling adi- 
abatic passage into a direct exciton. Depopulation of the 
output state occurs into the vacuum state. 



B. Qubit rotation 



In order to perform a desired rotation operation in- 
volving input states Ui-X") and \iqX), the system is adi- 
abatically driven into the central anticrossing occurring 
at Fr ~ 2.3kV/cm, see central box region in Fig. 3a. 
Figure [4k shows a sweep that drives the system from the 
anticrossing at Fj ~ 43.4kV/cm to the one at Fr. There, 



the system evolves by its internal dynamics. This is em- 
phasized by the absence of the coupling 51 in U, and the 
switch off of the laser once the qubit is initialized. The 
coherent oscillations allow for qubit rotations in the pi- 
cosecond scalepSwith a characteristic time ^j = 91.8ps, 
and exhibit a near unitary amplitude within the time 
frame shown in Fig. 4d. Notice that the coherent oscilla- 
tion relaxes on a longer time scale. This is due to strong 
direct exciton relaxation rates, and the weak effects of 
exciton virtual transitions occurring outside the qubit 
subspace, on the matrix element U{z) and energy shift 
of the indirect exciton, 5i(z), see Eqs. (12 1 and (J9J), re- 
spectively. The duration of this rotation determines how 
much population is transferred into the output state, in 
other words, how much the final indirect molecular state 
would follow the eigenvalue line, liJ-X"), upon bias rever- 
sal. 



C. Readout 

Once the qubit rotation has taken place, and the pop- 
ulation is transferred into the output state after a 7r rota- 
tion, the molecular eigenstate follows a different running 
eigenvalue in reverse bias. This is observed in Fig. 4e af- 
ter t = 3ns; here the applied bias pulse drives the output 
state lio-X") along the dressed spectral line \iqX), start- 
ing at Fj ~ 2.3 and finishing at F R ~ 43.4kV/cm. The 
fidelity of the read out de pends on a conditional adiabatic 
population passage^EIl from US-X") into the direct exci- 
ton lio-^O) then with the partial population transfer into 
the vacuum \qqX) (green solid and black dashed line, Fig. 
4e) . At the end of the sweep, far away from the central 
anticrossing, the direct exciton \\qX) is depopulated, by 
recombination emitting luminescence, without perturb- 
ing the indirect states or any other nearby exciton. 



D. Stability of coherent control 

We should emphasize that the control scheme requires 
the central indirect-exciton anticrossing to be isolated 
from other exciton states and resonances associated with 
transitions out of the qubit subspace. For a wide range 
of system parameters, the coherent rotation regime is 
achieved via an anticrossing which appears isolated in 
an energetically narrow window (easily detunable from 
other transitions), and occurs even if the ground states 
of the QDM dots are non-resonant ( S\q ^ Sq\ ) , protect- 
ing the qubit subspace and enhancing coherence. If other 
excited states (such as those associated with excited elec- 
tron/hole levels of the molecule) appear in the vicinity of 
the qubit window, their effects can be naturally incor- 
porated in the description. They may result in changes 
of the initialization field and pumping, but they would 
not intrinsically affect the main qubit rotation scheme. 
Certainly strong distortion of the relevant anticrossings, 
Figs. 3c-d, by a nearby state affects the rotation and ini- 



tialization fidelity since in that case the system would 
not be approximated by a well-separated two-level sys- 
tem. Further complications could arise if charge tunnel- 
ing rates into the diode contacts compete with the con- 
trol time scales; in that case, charged excitons (negative 
and positive trions) would not be negligible, affecting the 
charge stability of the exciton qubits, and enabling deco- 
herence by exciton spin dephasing mediated by electron- 
hole exchange interaction.^ However, these constraints 
can be relaxed by proper geometrical engineering of the 
excitonic spectrum and excitation conditions, and by se- 
lection of suitable molecules among the many produced 
in typical processes. 



VI. DISSIPATION EFFECTS 

We have assumed that the neutral indirect excitons are 
optically inactive for the chosen interdot distanced So 
that their intrinsic recombination rate is Ti = 0. How- 
ever, we have found that even with this assumption, the 
effect of the direct exciton spontaneous recombination, 
plus the influence of virtual transitions (mediated by tun- 
neling) occurring outside the qubit subspace, provides 
the molecular indirect exciton with a finite effective os- 
cillator strength and lifetime, which is ultimately a con- 
sequence of interdot quantum coupling. 
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FIG. 5. (Color online) Population relaxation dynamics for 
indirect excitons, from the numerical solutions of Eq. 4 with 
up to 14 exciton states included, (a) Population decay af- 
ter switching off the laser light in the initialization regime 
at Fi = 43.4kV/cm. Dashed line indicates the case where 
biexciton states are taken into account. Inset shows the cor- 
responding depopulation into the vacuum, (b) QDM internal 
dynamics in absence of optical perturbations, shows decay of 
Rabi oscillations of \o{X) and |?o-X} (red line, blue line) in the 
qubit rotation regime at Fa — 2.3kV/cm. Inset shows Rabi 
flops in the early stage of the dynamics 

Figure 5a shows the population time dependence of 



the input state \%X) after switching off the excitation 
power, £1(t) = 0, at a time r = 200ps, as in Fig. 4c 
(from the numerical solution of Eq. 4, with up to 14 ex- 
citon states included) . We see that the population relaxes 
into the vacuum (see inset) with a lifetime TJ ~ 218ns. 
Comparatively the dashed line indicates depopulation of 
the state when the biexcitonic degrees of freedom (and 
all 14 exciton basis states) are taken into account (no- 
tice they have no influence on the relaxation time when 
£1 = 0). Figure 5b shows the population time dependence 
of the QDM, when the input state has been driven adi- 
abatically into the rotation regime. Here the system is 
driven solely by its internal dynamics, which arises due to 
strong electron and hole tunneling, and the small influ- 
ence of FRET. Both indirect excitons (blue and red solid 
lines) enter a coherent oscillation regime, relaxing into 
the vacuum with a lifetime TJ ~ 25.2ns, with vanish- 
ing population transfer into the direct excitons. Notice 
that the relevant coherent oscillation period (~ 91.8ps), 
is orders of magnitude shorter than the relaxation time. 
On the other hand, our results indicate that despite as- 
suming infinite lifetime for the bare indirect excitons at 
very large bias, |F| ^ lkV/cm, the interdot couplings 
provide a strong dissipation channel for indirect excitons 
to relax into the vacuum at any finite value of F. In fact, 
experimental results have shown electric field tuning of 
radiative lifetimes at the direct-indirect molecular exci- 
ton resonance (near the tunneling induced anticrossing, 
F ~ ±20kV/cm), still in the range of ~ 2ns to ~ 10ns, if 
acoustic phonon mediated inter-level relaxation processes 
are present P3 

Calculation of the population damping of the state 
loi-X") inside the effective subspace Fa, as well as the 
damping of coherent oscillations within the qubit sub- 
space Vi, gives the effective lifetime of the qubit input 
state during initialization and the decoherence time dur- 
ing qubit rotation, respectively. Although we consider 
just the effects of intrinsic spontaneous recombination of 
the direct excitons, the effective decay rates of the molec- 
ular states can be significatively different, since they de- 
pend strongly on all interdot coupling mechanisms, and 
thus can be tuned at will as function of laser detuning, 
pump power and applied electric field. 

In what follows we are interested in the dissipative dy- 
namics of the indirect excitonic subspace. Since Vp <C 
|A|, we can ignore the effects of FRET in the following 
discussion. In order to obtain an analytical expression for 

the effective decay rates, we start with the non-Hermitian 
HamiltoniarPml 



which must of course reduce to Eq. ( 13 1 for Tx = Vp = 



H r = H- 



iT 



x n w 



(\%X)(™X\ + \%X)(%X\), (14) 



and project it onto the subspaces Va and Vp Conse- 
quently the resultant effective Hamiltonian, iFj {z), will 
be non-Hermitian, with matrix elements that depend on 
the direct exciton decay, Tx- Projection onto the sub- 
space Vp results in 

H?\z) = H<p + H£(z) + iH { /2(z) , (15) 



The diagonal part of the non-Hermitian projected Hamil- 
tonian is, 



\ 01 

with a projected real perturbation given by 



m^)= 



r^+4(z-A) 2 \v{z) a(z) J ' 



a(z) v{z) 



(16) 



(17) 



and its corresponding imaginary perturbation given by 



1 



^ = r| + 4 ( ;-A)4 7f'- ,llS1 



P 7 



with matrix elements given by, 

a(z) = m + t 2 h )(z-A), /3 = -2(ti + tl)T 



x 



v(z) = 8t e t h (z-A), j = -4t e t h T 



(19) 



x 




FIG. 6. (Color online) Relaxation times for bare and dressed 
indirect excitons as function of applied electric field F. Light 
red (blue) line corresponds to relaxation time for the dressed 
(bare) indirect exciton eigenstates. Far from anticrossing 
value of electric field Fi = 2.38kV/cm, the dressed relaxation 
time approaches the bare exciton relaxation time as molecular 



indirect eigenstates of H r 
\11X)- 



(i) 



become the indirect states |ioX), 



The dissipative part of the Hamiltonian, -£/},„, con- 
tains two sources of decoherence, a diagonal source pro- 
portional to /3 and an off-diagonal source proportional 
to 7. P is the intrinsic population relaxation rate of the 
dressed indirect excitons, which is non-vanishing when ei- 
ther the electron or the hole tunneling is non-zero. This 
term is the dominant part of the relaxation rate, and is 
present along the whole electric field sweep, inside and 
outside the qubit subspace. 7 arises from the simultane- 
ous tunneling of an electron and a hole, and can be inter- 
preted as an interference term between the simultaneous 
paths of these two charges. This term is much smaller 



(as t e S> th), and contributes mostly inside the qubit 
subspace. Therefore, the damping of coherent oscilla- 
tions during the rotation operation, contains population 
decay and dephasing. 

The decay rates of the dressed states inside the qubit 
subspace, Vx, can be obtained by diagonalization of 
the Hamiltonian (15), which gives the total imaginary 
component of the indirect dressed excitons. This to- 



tal rate can be written as Vd(z) = Im[Diag(_ffp ')]. 
= r/(z) + r$(z), for either indirect exciton \j), where 



f/W 



m+ti)r 



x 



r 2 

1 X 



4(z-A) 2 



(20) 
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+ A, 
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2\2 
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(21) 
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(r 2 
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Figure 6 shows the dependence of the effective relax- 
ation times, FJ (z) (blue line) and T^ (z) (light red line), 
as function of applied electric field F. For values of field 
away from the central anticrossing in Fig. 3c, the effec- 
tive relaxation time, of the eigenstates of Hj, (z), TJ- l 1 (z), 
approaches the relaxation time of the bare effective in- 
direct excitons, TJ (z), where the interference term (de- 
phasing), f$, significatively diminishes. At the avoided 
crossing with Fr — 2.3kV/cm, the eigenstate relaxation 
time reaches a minimum value of ~ 22.2 ns. Notice, it is 
Tj (z) which is more physically relevant, as it gives the 
lifetime of the qubit logical states for all values of F, be- 
tween the initialization regime, at Fj = 43.4kV/cm, and 
the rotation regime at Fr = 2.3kV/cm, with correspond- 
ing values of 214ns and 24.5ns, in good agreement with 
the numerical solutions of the Lindblad equation shown 
on Fig. 4. Notice that the ratio between gate rotation 
and decoherence times, U/Tj, is monotonically depen- 
dent on the interdot distance, so that an optimal ratio 
is not given within the model, but is set by the QDM 
geometry. 



VII. CONCLUDING REMARKS 



r 



ton relaxation rates arising from spontaneous decay, the 
subspace of indirect qubit states has large decoherence 
times. This is explained by a large suppression of the in- 
teraction between the qubit subspace and the laser field. 
On the other hand, the qubit can be initialized with near 
unity fidelity via higher order couplings to the radiation 
field. In this manner the input state can be shelved for 
time intervals well beyond the direct exciton relaxation 
time. The use of an adiabatic bias pulse permits driving 
the input state into different molecular resonances, in 
particular, into a resonance that mixes coherently the in- 
put state with a target state of the qubit. Interestingly a 
reverse bias pulse drives an adiabatic passage of the out- 
put logical state through a tunneling induced anticross- 
ing, transferring half of its population to a spatially direct 
exciton. This enables the possibility of directly reading 
the output of the qubit rotation. The suppression and 
tunability of the exciton interactions with states outside 
the qubit subspace contrast drastically with qubits de- 
fined via excitons with spatially direct character. This 
opens the possibility of using neutral exciton states as 
elemental blocks within a complex QDM quantum com- 
putation scheme that uses indirect excitons, such as spin 
qubits in molecular trions. 



In summary, we have shown that the exciton spectrum 
of a QDM can be used to define an optimally defined 
qubit using two spatially indirect neutral excitons. We 
found that the interplay of optical excitation and charge 
tunneling can produce optical signatures that identify the 
indirect exciton qubit subspace. Although the QDM is 
treated explicitly as an open quantum system, with exci- 
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